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Objectives 

This  research  program  is  focused  on  the  study  of  the  atomic  force  microscope  (AFM)  as 
an  instrument  for  quantitative  measuring  and  imaging  material  properties  with  nanoscale 
resolution.  More  specific  to  this  research  are  the  dynamic  modes  of  AFM,  in  which  the 
specimen  moves  relative  to  the  AFM  cantilever  tip  while  in  contact.  Particular  emphases 
include  the  higher-order  vibration  modes  and  the  nonlinear  vibration  response.  In 
addition  to  the  original  objective,  improved  modeling  of  AFM  .cantilevers  has  been 
undertaken.  The  use  of  simplified  beam  models  with  uniform  cross  section  does  not 
reflect  the  AFM  cantilevers  that  are  often  used  in  experiments.  High  resolution  cantilever 
measurements  with  a  scanning  electron  microscope  (SEM)  and  modeling  with  finite 
element  analysis  have  been  used  to  improve  the  comparisons  between  experiments  and 
numerical  models. 


Description  of  Effort 

During  this  project,  the  development  of  AFM  experimental  expertise  was  undertaken.  An 
AFM  was  purchased  and  modified  for  the  dynamic  AFM  measurements.  The  system 
allows  the  cantilever  resonance  frequencies  to  be  determined  both  in  and  out  of  contact 
with  a  specimen.  In  addition,  the  dynamic  implementation  allows  the  surface  topography 
and  material  stiffness  to  be  mapped  simultaneously.  During  this  project,  several  aspects 
of  this  research  were  explored.  One  focus  was  on  comparison  of  AFM  methods  with 
other  experimental  measurement  techniques  using  thin  films  as  a  test  base.  Improved 
modeling  of  the  AFM  cantilevers  was  undertaken  as  it  became  apparent  that  simpler 
methods  were  not  sufficient  for  experimental  analysis.  Finite  element  models  are  being 
constructed  to  include  observed  variations  in  cantilever  width,  thickness,  and  tip  location. 
The  analytical  modeling  has  been  focused  on  the  nonlinear  vibrations  of  AFM  cantilevers 
with  nonlinear  boundary  conditions.  The  studies  of  nonlinear  vibrations  have  been 
extended  for  contact  boundary  conditions  of  a  general  nature  such  that  both  soft  and  hard 
contacts  may  be  analyzed.  The  method  of  multiple  scales  is  the  primary  method  used  for 
these  studies.  Finally,  the  coupling  of  flexural  and  torsional  vibration  modes  has  been 
examined.  The  coupling  is  the  result  of  the  specific  geometry  of  certain  AFM  cantilevers. 
Both  experimental  and  analytical  studies  have  been  conducted. 
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Accomplishments 


Comparisons  with  other  Techniques.  The  experimental  setup  for  the  dynamic  AFM 
experiments  (also  known  as  atomic  force  acoustic  microscopy  or  AFAM)  is  shown 
schematically  in  Fig.  1.  A  commercial  AFM  (Thermomicroscopes  Autoprobe  CP)  was 
modified  for  these  experiments.  A  reference  specimen  and  an  unknown  specimen  are 
placed  on  a  broadband  ultrasonic  transducer  (2.25  MHz  center  frequency).  The 
transducer  is  mounted  on  the  AFM  stage  and  is  excited  by  the  function  generator.  The 
output  from  the  AFM 
photodiode  is  mixed 
with  the  driving  signal 
in  the  lock-in  amplifier. 

The  lock-in  output  and 
the  photodiode  output 
are  both  used  by  the 
AFM  for  imaging.  The 
function  generator  is 
controlled  by  Labview 
software. 


Figure  1.  Schematic  of  the  atomic  force  acoustic  microscopy 
experiment.  Both  the  topography  and  lock-in  amplifier  output  are 
used  for  imaging. 


In  a  typical  experiment,  the  free  and  contact  cantilever  resonances  are  measured 
on  both  the  reference  and  the  sample  of  interest.  The  cantilever  resonances  are  found  by 
scanning  an  appropriate  frequency  range  with  the  function  generator  and  plotting  the 
lock-in  output.  The  contact  resonances  are  functions  of  applied  load  on  the  AFM 
cantilever.  The  frequency  shifts  may  be  used  to  determine  the  mechanical  properties  of 
the  unknown  relative  to  the  reference.  First,  the  necessary  spring  stiffness  of  the  contact 
k  is  determined  from  the  natural  frequency  data  and  an  appropriate  vibration  model.  The 
value  of  k  is  then  used  to  determine  the  contact  modulus  of  the  sample  using  the  Hertzian 

contact  relation  E*  =  E*(k*  / k*J'2 ,  where  the  subscripts  denote  the  sample  and  reference. 

The  sample  indentation  modulus  Ms  is  then  determined  from  (is’jr1  =  (m,(>  )"'  +  (Ms )"' , 

where  M,ip  is  the  indentation  modulus  of  the  tip.  AFAM  results  for  indentation  modulus 
on  samples  of  aluminum  and  niobium  were  compared  with  results  using  surface  acoustic 
waves  (SAW)  and  instrumented  indentation  (HT)  [1],  The  results  are  summarized  in 
Table  1.  These  results  are  typical  of  other  measurements  on  different  materials.  The 
primary  cause  of  the  wide  range  in  the  AFAM  results  is  thought  to  come  from  the  choice 
of  the  reference  material.  For 
the  results  in  Table  1,  both 
glass  and  single  crystal 
silicon  were  used  as 
references.  The  appropriate 
use  of  reference  samples  for 
AFAM  is  the  subject  of  ongoing  research. 


Table  1.  Comparison  of  indentation  modulus  M  (GPa)  for 
aluminum  and  niobium  using  several  different  methods  [11. 


Material 

Literature 

SAW 

IIT 

AFAM 

A1 

76-81 

78  ±1 

86  ±4 

54-81 

Nb 

116-133 

120  ±5 

— 

91  - 137 
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Imaging.  Once  the  resonances  are  measured,  an  AFAM  image  may  be  made.  The 
transducer  motion  is  fixed  at  a  single  frequency  slightly  higher  or  lower  than  the 
resonance  for  a  given  applied  load.  The  feedback  control  loop  of  the  AFM  is  set  to 
constant  force  and  the  scan  is  made.  The  topography  and  stiffness  images  are  obtained 
simultaneously.  Example  results  for  polysilicon  are  shown  in  Fig.  2.  The  topography 


Topography,  0518T015.HDF 


AU X4.  05180014. HDF 


Figure  2.  Example  results  of  atomic  force  acoustic  microscopy  (AFAM)  on  a  sample  of  polysilicon 
MEMS.  The  two  lefthand  images  are  topography,  while  the  two  righthand  images  are  the  AFAM 
images.  Of  particular  note  are  the  stiffness  contrasts  observed  within  single  grains. 

images  are  in  the  lefthand  column  while  the  stiffness  images  (lock-in  output)  are  shown 
in  the  righthand  column.  Results  from  two  scans  are  shown.  The  stiffness  images  for 
these  samples  show  the  interesting  contrast  within  regions  that  appear  to  be  single  grains 
in  the  topography  image.  Also  of  note  are  the  single  grains  that  appear  to  have  different 
stiffness.  Calibrated  stiffness  measurements  may  be  made  by  comparing  spectra 
measured  from  a  known  material,  such  as  single  crystal  silicon,  to  spectra  from  the 
unknown. 


Improved  Modeling.  Improved  modeling  of  AFM  cantilevers  was  undertaken  as  it 
became  clear  that  the  simple  models  based  on  cantilevers  of  uniform  cross-section  did  not 
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apply  to  many  commercially  available  cantilevers.  To  date,  AFAM  experiments  have 
been  interpreted  exclusively  with  an  analytical  model  in  which  it  is  assumed  that  the 
cantilever  cross-section  is  uniform  over  the  length  of  the  beam.  The  assumption  is  not 
strictly  true,  even  for  cantilevers  that  are  uniform  over  the  majority  of  their  length.  Thus, 
data  interpretation  using  such  analytical  models  is  an  approximation  for  real  experimental 
conditions.  Such  approximations  are  very  good  in  some  cases.  However,  cantilevers  such 
as  dagger  shaped,  V-shaped,  or  triangular  cantilevers  are  useful  for  measurements  but  are 
clearly  not  uniform. 

To  address  the  current  deficiency  in  these  simplified  AFAM  models,  finite 
element  models  (FEM)  have  been  developed  for  AFAM  applications.  Variations  in 
geometry  that  extend  over  a  small  range  of  the  cantilever,  such  as  the  triangular  portion 
of  dagger-shaped  cantilevers,  are  easily  modeled  using  FEM.  Two  types  of  cantilevers 
were  investigated,  one  with  rectangular  plan  view  and  one  with  a  dagger-shaped  plan 
view  [1],  For  both  cantilevers,  the  length  and  width  dimensions  were  measured  using  an 
optical  microscope.  The  thickness  was  measured  using  scanning  electron  microscopy 
(SEM).  The  optical  measurements  were  more  accurate  (relative  to  the  degree  of  accuracy 
required)  than  the  SEM  measurements  due  to  mounting  and  alignment  issues.  Moreover, 
the  beam  resonance  frequencies  are  much  more  sensitive  to  thickness  than  to  length  or 
width.  For  these  reasons,  the  cantilever  thickness  was  the  primary  variable  adjusted  to 
match  the  experimental  results.  Other  variations  in  the  different  ^cantilevers  were  also 
included  in  the  models.  For  instance,  the  tip  of  the  rectangular  cantilever  was  not  located 
exactly  at  its  end,  but  some  unknown  distance.  Therefore,  the  tip  position  was  allowed  to 
vary  in  the  FEM.  The  location  of  the  dagger  cantilever’s  tip,  on  the  other  hand,  was 
exactly  at  the  end  and  was  not  a  variable  parameter  in  the  model.  Other  parameters  that 
were  investigated  in  the  FEM  include  tip  mass,  damping,  and  cantilever  thickness 
variations.  Example  numerical  results  for  both  types  of  cantilevers  are  given  in  Table  2. 
The  indentation  modulus  was  calculated  from  the  contact  frequencies  for  both  aluminum 
and  niobium.  The  results 
are  based  on  two  reference 
materials.  The  FEM  results 
in  some  cases  are  not 
necessarily  an  improvement 
over  the  simplified  model. 

However,  it  should  be 
noted  that  in  many  cases 
unphysical  parameters  are 
used  in  the  simplified 
models,  such  as  tip 
placement  very  far  from  the 
beam  end.  Further 
improvements  to  these 
modeling  efforts  are  being 
examined. 


Table  2.  Comparison  of  results  for  indentation  modulus  M  using 
analytical  and  finite  element  models  of  AFM  beams  jl]. 


Material 

|  Cantilever  |  Reference 

Model 

M 

A1 

rectangular 

glass 

uniform 

76  ±3 

FEM 

76  ±10 

A1 

dagger 

glass 

uniform 

64  ±3 

FEM 

54  ±4 

Nb 

rectangular 

glass 

uniform 

89  ±7 

Si 

140  ±6 

glass 

FEM 

95  ±  14 

Si 

128  ±  16 

Nb 

dagger 

glass 

uniform 

87  ±  1 

Si 

139  ±  5 

glass 

FEM 

92  ±  6 

Si 

120  ±14 
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Nonlinear  Vibrations.  The  nonlinear  forces  that  develop  between  the  AFM  tip  and  the 
specimen  often  result  in  nonlinear  behavior  in  the  vibrations.  Experimental  results  have 
shown  the  type  of  softening  nonlinearity  expected.  The  fundamental  problem  of  the 
nonlinear  vibrations  of  a  beam  with  cantilever-Hertzian  contact  boundary  conditions  was 
investigated  [2].  The  method  of  multiple  scales  was  used  to  analyze  this  problem  in 
which  it  was  assumed  that  the  beam  remains  in  contact  with  the  moving  surface  at  all 
times.  The  contact  between  the  AFM  tip  and  the  sample  surface  was  modeled  as 
Hertzian.  One  primary  result  from  this  analysis  is  the  amplitude-frequency  relation  for 
the  various  flexural  modes.  The  amplitude-frequency  curves  exhibit  softening  behavior  as 
expected.  The  amount  of  softening  depends  on  the  linear  contact  stiffness  relative  to  the 
beam  stiffness  (characterized  by  dimensionless  parameter  /?)  as  well  as  the  specific  mode. 
An  example  result  is  shown  in  Fig.  3.  The  modal  amplitude  is  plotted  versus  the  detuning 
parameter  for  the  first  and  second  modes  for  two  values  of  p.  The  sensitivity  to  the 
nonlinearity  is  similar  to  the  linear  sensitivity  [3].  The  sensitivity  of  the  different  modes 
is  reflected  in  their 
nonlinear  behavior  — 
some  modes  exhibit 
larger  shifts  in  natural 
frequency  than  others 
for  the  same  type  of  | 
contact  stiffness.  The  t 
associated  nonlinear 
normal  modes  of  this 
system  were  also 

derived.  The  modes 

include  a  nonlinear 

addition  to  the  linear,  Figure  3.  Nonlinear  frequency  shift  of  the  first  and  second  modes  for 

harmonic  component  as  two  values  of  linear  contact  stiffness  ft  =  50  (solid)  and  500  (dashed), 

well  as  a  static  offset  Each  mode  responds  differently  to  the  nonlinearity  [2j. 

term  and  second  and  third  harmonic  components. 


Detuning  Parameter,  a  Detuning  Parameter,  <r 


Modal  Coupling.  The  coupling  of  flexural  and  torsional  vibration  modes  of  atomic  force 
microscope  cantilevers  has  also  been  investigated.  This  type  of  coupling  occurs  in 
cantilevers  of  certain  geometries.  Solutions  were  obtained  using  an  analytical  bending 
and  torsional  model  as  well  as  using  the  finite  element  method.  The  coupling  considered 
arises  from  small  geometric  asymmetries  in  the  cantilever,  such  as  an  offset  tip  or  uneven 
distribution  of  cantilever  mass.  The  asymmetries  lead  to  coupling  of  the  respective 
flexural  and  torsional  governing  equations.  Specific  conditions  necessary  for  the  coupling 
of  the  first  torsional  mode  with  various  flexural  modes  were  derived.  For  isotropic 
rectangular  cantilevers  with  rectangular  cross-section,  the  tendency  for  coupling  is 


characterized  by  a  dimensionless  parameter  H = 


,  where  w  is  the  width,  L  is 


the  length,  and  v  is  Poisson’s  ratio.  Conditions  for  coupling  are  easily  satisfied  by 
commercially  available  AFM  cantilevers.  Example  experimental  results  are  shown  in  Fig. 
4  [4].  The  experiments  were  done  with  commercial  rectangular  AFM  cantilever  beams 
made  of  silicon.  The  cantilevers  were  set  into  vibration  using  a  piezoelectric  ultrasonic 
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transducer  coupled  to  the  chip  of  a  cantilever.  The  transducer  was  excited  harmonically  in 
the  frequency  range  of  100  kHz  to  3  MHz.  The  vibrations  were  monitored  using  an 
optical  Michelson  heterodyne-interferometer  allowing  the  surface  of  the  cantilever  under 
examination  to  be  scanned  with  a  lateral  resolution  of  several  microns.  The  measured 
frequency  spectrum  was  used  to  identify  the  various  modes.  The  amplitude  and  phase  of 
four  free  vibration  modes  of  the  beams  are  shown  in  Fig.  4.  Theoretical  predictions  of  the 
modal  coupling  are  shown  in  Fig.  5  for  H=  0.025.  These  results  agree  very  well  with  the 
experimental  results. 


Exciting  frequency  21 1 .81  kHz,  Uss=  6  V 


Deflection 

NO  nm 
9  nm 
18  nm 


Exciting  frequency  26539  kHz,  Uss=  6  V 


BO  nm 
16  nm 
31  nm 


Fxritinn  frpnupnrv  Q  kW?  (i  =R\/ 


-  n 
0 
71 


Figure  4.  Amplitude  and  phase  distribution  of  the  first  four  modes  of  a  rectangular  Si-cantilever  that 
are  predominantly  torsional.  The  scan  size  is  about  550  x  100  pm2  [4]. 


Figure  5.  Predicted  third  and  fourth  mode  shapes  calculated  from  the 
flexural/torsional  coupling  model  with  JJ  =  0.025.  The  interaction  observed 
is  a  coupling  between  the  first  torsional  and  third  flexural  modes  [4]. 


Interactions.  Collaborative  efforts  have  continued  with  Dr.  D.  C  Hurley  of  the  Materials 
Reliability  Division  of  the  National  Institute  for  Standards  and  Technology  (NIST)  in 
Boulder,  Colorado.  This  collaboration  includes  the  sharing  of  test  samples  and  AFM 
cantilevers,  comparison  of  experimental  data,  and  collaborative  theoretical/experimental 
research.  Collaborative  research  is  also  continuing  with  the  group  headed  by  Dr.  W. 
Arnold  at  the  Fraunhofer  Institut  fur  zerstorungsfreie  Prufverfarhen  (IZFP)  in 
Saarbriicken,  Germany.  The  PI  and  students  are  collaborating  with  the  IZFP  on  the 
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development  of  finite  element  models,  analytical  vibration  models,  and  experimental 
comparisons.  Collaboration  with  the  thin  films  group  at  UNL  (S.  Rohde)  has  also  proven 
fruitful  [5].  The  AFAM  research  is  continued  through  support  of  the  National  Science 
Foundation  (ceramic  nanofibers)  and  the  Department  of  Energy  (sintered  materials). 
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Transitions 


As  an  ongoing  part  of  this  research,  the  analysis  and  calculations  used  for  AFAM  have 
been  developed  into  web-based  calculation  tools  ('http://em-iatumer.unl.edu/calcs.htm V 
For  example,  one  of  these  tools  calculates  AFM  cantilever  frequency  spectra  and  modal 
sensitivity  as  a  function  of  the  surface  stiffness  for  both  flexural  and  torsional  vibrations 
(shown  in  Fig.  6).  Current  options  include  three  cantilever  plan  views  of  rectangular,  V- 
shaped,  and  triangular.  The  user  chooses  the  cantilever  shape  and  inputs  material  property 
and  geometric  data.  The  information  is  then  passed  to  a  MATLAB  program  that  performs 
the  calculations.  The  output  is  then  posted  in  plot  or  table  format  in  the  browser  frame. 
Additional  tools  focused  on  AFAM  experiments  are  also  now  included. 


Figure  6.  Example  output  from  one  of  the  web-based  AFM  calculation  took. 
Several  different  cantilever  geometries  are  included  as  are  flexural  and 

torsional  modes. 
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Commercial  rectangular  atomic  force  microscope  cantilever  beams  made  of  silicon  were  set  into  vibration,  using  a 
piezoelectric  ultrasonic  transducer  coupled  to  the  chip  of  a  cantilever.  The  transducer  was  excited  with  continuous  rf  in 
the  frequency  range  of  100  kHz  to  3  MHz.  The  vibrations  were  monitored  using  an  optical  Michelson  heterodyne- 
interferometer  allowing  the  surface  of  the  cantilever  under  examination  to  be  scanned  with  a  lateral  resolution  of  several 
pm.  A  number  of  free  torsional  and  flexural  vibration  modes the  beams  were  imaged  quantitatively.  Comparison  of 
the  experimental  resonance  frequencies  and  the  amplitude  and  phase  distribution  of  the  modes  to  theoretical  models 
showed  that  asymmetries  in  the  beam  strongly  influence  the  vibrational  behavior  of  the  beam.  The  consequences  for 
quantitative  local  stiffness  measurements  are- discusSbd.: : 
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imaging  of  vibration  a  F 


22  1.  Introduction 


Since  the  invention  bf  the  Atomic  force  micro¬ 
scope  (AFM)  [1]  several  different  quasi-static  and 
dynamic  operating  modes  have  been  developed, 
for  example  frietioil  force  microscopy  (FFM)  [2] 
or  atomic  force  acoustrc  imcroscopy  [3].  During 
the  past  few  yeais*  rpany  groups  have  been  im- 
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68 1  -93025903^^1? :  •  v  j 

E-mail  address:  reinstaedtlemi2fp.fhg.de  (M.  Reinstacd- 
tler). 


proving  the  quasi-static  FFM  technique  [4,5]  to  a  29 
dynamic  method,  expanding  the  possibility  to  ex-  30 
amine  the  origin  of  friction  on  a  nanoscale  [6-18]  31 

without  the  shortcomings  of  conventional  FFM,  32 
which  are  time-consumptive  measurements  with  33 
high  noise  and  unstable  friction  imaging.  34 

One  crucial  component  for  all  AFM -operating  35 
methods  is  the  flexible  force-sensing  cantilever  [19].  36 

These  micro-fabricated  AFM  cantilevers  can  be  37 
regarded  as  miniature  elastic  beams  with  charac-  38 
teristic  bending  and  torsional  vibration  modes.  By  39 
means  of  elasticity  theory  it  is  possible  to  calculate  40 
the  resonance  frequencies  and  the  local  vibration  41 
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42  amplitudes  and  phases  of  the  modes.  Depending 

43  on  the  geometry  of  the  beam  analytical  or  nu- 

44  merical  methods  [20,21]  can  be  used.  However,  the 

45  exact  cantilever  geometry  is  often  not  known  and 

46  can  vary  considerably  within  one  set  of  comraer- 

47  cial  beams.  Properties  of  the  beams,  such  as  ver- 

48  tical  and  lateral  spring  constants  are  of 

49  fundamental  importance  for  quantitative  mea- 

50  surements,  in  particular  for  local  elasticity  mea- 

51  surements  [20].  By  examining  the  vibration  spectra 

52  of  the  cantilevers  during  our  ultrasonic  AFM 

53  measurements,  we  often  noticed  the  appearance  of 

54  additional  resonances  not  expected  theoretically. 

55  Here  we  show  that  the  free  vibrational  modes  of 

56  clamped- free  silicon  AFM  cantilevers  can  be  im- 

57  aged  quantitatively.  This  provides  evidence  about 

58  the  origin  of  the  additional  modes  and  also  allows 

59  a  quality  control  of  AFM  sensors  because  the  ef- 

60  fective  spring  constant  depends  on  the  vibrational 

61  mode  of  the  cantilever  [22],  which  in  turn  is  in- 

62  fluenced  by  the  geometry  of  the  cantilever,  the 

63  mass  distribution  along  its  length,  and  the  position 

64  of  the  tip. 


65  2.  Theory 

66  2.J.  Flexural  vibrations 


67 

68 

69 

70 

71 

72 
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78 
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81 


Bending  and  torsional  vibrations  of  uniform, 
homogenous  beams  with  constant  ^rectangular 
cross-section  can  be  described  by  a  four^-  and  a 
second-order  partial  differential  equation,  respec¬ 
tively.  The  equation  of  motion  for  the' bending 
modes  can  be  found  in  textbooks  [23] 


04y  02y 

£/_+p/,_  =  0 


(i) 


where  E  is  the  Young's  modulus-  7  =  abl  j\2  is  the 
area  moment  of  inertia  of  a  beam  with  rectangular 
cross  section,  a  is  itawidth  and  b  its  thickness,  p  is 
the  mass  density  and  A  =  ab  is  the  cross  sectional 
area.  Here,  x  "^Ttne  ^iordinate  in  the  longitudinal 
direction  of  the  cantilever,  j>(x)  is  the  excursion 
from  the  rest  position  of  the  length  element  at  x.  A 
general  solution  -for  Eq.  (1)  can  be  written  as 

y(x,t)  =  («,e^+  a2e->*  +  a^x  +  me-*)#*  (2) 


Here,  co  =  2nf  is  the  angular  frequency  and 
y  =  2k/ X  the  flexural  wave  number.  Inserting  Eq. 
(2)  in  Eq.  (1)  yields  the  dispersion  relation 


-  pAar  =  0 


Ely* 

or 


C0  =  yV^A 


y;:. 


(3) 


(4) 


The  constants  a,  (/  =  determined 

by  the  boundary  conditions^  free  vibrations,  the 
cantilever  is  clamped  atgc;  —  0,  \yhereas  the  free  end 
is  at  x  —  L.  Thiryidlds  Me  Characteristic  equation 
for  the  wave  numberi^  .  -;.:; 

cosynLcoshy^^  1  *  (5) 

Using  Eq.  (4)%id  the  solutions  y „L  of  Eq.  (5),  one 
obtains  the  bendmg-mode  eigenfrequencies 


0% 


(6) 


a  The  amplitude  distribution  of  the  flexural  can¬ 
tilever  vibration  y„  (x)  of  each  mode  is  described  by 
the  following  equation 


M*)  —yo  ' 


(cos  ynx  -  cosh  y,pc) 


cos  y„L  +  cosh  y„L 
sin  ynL  +  sinh  y„L 


(sin  yltx  +  sinh  y„x) 


(7) 


where  y^  is  the  vibrational  amplitude. 
2.2.  Torsional  vibrations 


The  equation  of  motion  for  torsional  modes  is 


02<9  rd2<9  A 

CTw-pJw-° 


(8) 


where  <9  is  the  angle  of  twist,  cT  =  ablG/ 3  is  the 
torsional  stiffness  and  J  =  a1  bj  12  the  polar  area 
moment  of  inertia.  A  general  solution  for  Eq.  (7) 
can  be  written  as 


0(x)  =  A  sin  t]x  -f  B  cos  rjx  (9) 

where  constants  A  and  B  are  also  given  by  the 
clamped-free  boundary  conditions  for  the  cantile¬ 
ver  oscillations,  so  the  amplitude  distribution  is 
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115 

116 


®(x)  —  A  sin  r\x 


(10) 


The  torsional  wave  number  7]  is  obtained  by  the 
balance  of  the  torsional  moments  [12,24] 


Cl 


d0(x) 

dx 


=  0 


118  The  eigenfrequencies  for  torsional  vibrations  are 

119  obtained  by  solving  Eq.  (8) 


2/i-l  1  b  Ig 

2  L  a]j~p 
ay 


(11) 


121  Unlike  the  bending  mode  eigenfrequencies  (Eq. 

122  (6))  the  torsional  modes  are  equally  spaced  by  A/, 

123  because  Eq.  (8)  exhibits  no  dispersion. 


124  3.  Experiments 

125  3.1.  Cantilever 

12 6  The  examined  cantilevers  were  manufactured 

127  from  highly  doped,  single-crystal  silicon  with  in- 

128  tegrated  single-crystal  silicon  tip.  The  tip  jls 

129  pointing  into  the  (100)  direction,  whereas  the 

130  cantilever  pointed  into  the  (110)  direction^  We 

131  used  two  low  force  constant  types  d^i^ed  Tor 

132  contact  mode  (Fig.  1)  and  for  lateraii  forde  mi- 

133  croscopy,  respectively,  from  the  same  riiairiufaC- 

134  turer  (Pointprobe®,  Nanosensorsil  Wetzlar, 

135  Germany).  For  generating  a  high  enough  torsional 

136  vibration  amplitude,  we  selected ; cantilevers  with 

137  low  bending  and  hence  also  low  torsionaf  stiffness. 

138  In  order  to  determine  the  geometrical  parameters, 

139  the  width  and  length  of  the  cantilevers  were  mea- 


Fig.  1 .  Optical  micrograph  of  the  cantilever  used  for  torsional 
mode  measurement  (type  CONT);  Calculated  spring  constant 
k  —  0.131  N/m;  nominal  thickness:  2  ±0.5  pm.  The  cantilever  is 
supported  by  a  single  crystal  silicon  holder  as  can  be  seen  on  the 
right. 


sured  with  an  optical  microscope.  From  the  mea¬ 
sured  resonance  frequency  of  the  first  ten  bending 
modes  the  thickness  of  the  cantilever  was  calcu¬ 
lated  using  Eq.  (6)  avoiding  time-consuming 
measurements  using  an  electron-microscope. 

3.2.  Atomic  force  microscope  set-up 

The  experimental  set-up  to  measure  the  free  vi¬ 
bration  spectra  of  the  cantilever  with  the  AFM  is 
shown  in  Fig.  2.  A  commercial  Scanning  Probe 
Microscope  (Dimension  3000,  Digital  Instru¬ 
ments,  Santa  Barbara,  :CA)  with  additional  ultra¬ 
sonic  equipment  was  used  to  excite  and  measure 
flexural  and  torsional  vibration  resonances  of  the 
cantilever.  A  frequency  generator  (HP  33120  A) 
provides  a  sinusoidal  signal,  which  is  applied  to  a 
conventional  ultrasonic  shear  wave  transducer 
(Panametrics  VI 51,  5  MHz  center  frequency)  at¬ 
tached  to  the  sample.  As  a  sample,  we  used  a  piece 
of  a  silicon  wafer.  The  ensuing  sample  surface  vi¬ 
bration  couples  via  the  air  column  between  the 
cantilever  and  the  sample-  into  the  cantilever 
mounted  over  the  sample  with  a  distance  of  about 
2-5  pm.  The  cantilever  vibrations  were  measured 
With  the  internal  optical  detector  of  the  instru¬ 
ment,  provided  its  bandwidth  was  in  the  MHz 
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Fig.  2.  Experimental  set-up  for  torsional  vibration  measure¬ 
ments  in  an  AFM. 


140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 


M.  Reinsiaedtler  et  al.  /  Surface  Science  xxx  (2003)  xxx-xxx 


regime,  otherwise  a  fast  external  amplifying  system 
was  used.  The  HF  amplitude  detection  scheme 
consists  of  a  heterodyne  down-converter  which  in 
addition  receives  a  reference  signal  at  the  excita¬ 
tion  frequency  from  the  frequency  generator.  The 
heterodyne  down-converter  operates  at  frequen¬ 
cies  from  75  kHz  to  10  MHz  and  shifts  the  desired 
signal  to  a  fixed  20  kHz  intermediate  frequency.  It 
is  evaluated  by  a  lock-in  amplifier.  The  set-up  is 
controlled  by  a  Labview  program  enabling  one  to 
change  the  frequency,  to  read  the  lock-in  output, 
and  to  store  the  spectra. 

Typical  measured  lateral  vibration  spectra  of 
cantilever  #1  (type  CONT)  are  shown  in  Fig.  3. 
Besides  torsional  modes  flexural  modes  also  ap¬ 
peared  in  the  lateral  channel,  however,  never  were 
torsional  modes  observed  in  the  vertical  channel. 
By  calculating  the  first  torsional  eigenfrequency 
using  Eq.  (11),  the  212  kHz  peak  in  the  lateral 
channel  spectrum  can  be  identified  as  belonging  to 
a  torsional  resonance.  However,  a  second  peak 
appeared  at  265.5  kHz  without  a  corresponding 
peak  in  the  vertical  spectrum  at  the  same  fre-  ; 
quency.  In  order  to  clarify  the  situation,  we  pe&: 
formed  interferometric  measurements  of  the"1" 
displacements  of  the  vibrating  cantilevers  in  .order 
to  compare  the  experimentally  determined  vibra¬ 
tion  patterns  to  the  expected  patterns  as  calculated 
by  Eqs.  (7)  and  (11).  ^  ♦ 

3.3.  Interferometric  measuring  system.  , 


The  interferometrical  set-up  is  shown  in  Fig.  4. 
The  chip  of  the  cantilever  Was  attached  with  a 
couplant  to  a  shear  ultrasonic  transducer  operated 
between  100  kHz  and  3  MHz.  It  emitted  into  the 
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Fig.  4.  Experimental  set-up  for  imaging  the  resonance  modes 
with  an  optical  interferometer. 

chip  transverse  .waves  polarized  perpendicular  to 
the  long  axis  of  thecahtileyer  beam  and  parallel  to 
its  width  causing  forced  •  vibrations.  The  absolute 
amplitudes  of  the  cantilever  vibrations  were  mea¬ 
sured  with  an'  optical  Michelson  heterodyne  in¬ 
terferometer  [25]’ IThe  He-Ne  laser  beam  of  the 
interferometer  was  focused  onto  the  cantilever 
with  a  spot  diameter  of  2-5  pm.  Using  a  motor- 
driven  linear  translation  stage,  the  focal  spot  of  the 
interferometer  was  scanned  over  the  surface  of  the 
cSbdleyir.  A  step  size  of  2  pm  was  chosen.  To 
determine  the  resonance  frequencies,  a  spectrum 
•::*yas  taken  at  a  location  close  to  the  free  end  of  the 
cantilever.  Then  the  modes  were  imaged  by  scan¬ 
ning  the  vibrating  cantilever  at  its  resonance  fre¬ 
quencies  and  detecting  the  deflection  of  each 
position  as  a  function  of  time.  The  amplitude  and 
phase  of  each  time  signal  were  obtained  by  using 
Digital  Fourier  Transformation. 

The  measured  amplitude  and  phase  distribution 
of  the  third  bending  mode  of  a  225  pm  long  beam 
are  shown  in  Fig.  5.  Two  nodes  and  three  crests  of 
the  vibration  are  found  according  to  theory.  The 
phase  of  the  signal  is  measured  relative  to  the  ex- 
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Fig.  3.  AFM  lateral  channel  spectra  showing  additional  modes. 
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240  cantilever  base.  This  is  caused  by  the  lower  bend- 

241  ing  stiffness  of  the  cantilever  compared  to  its  tor- 

242  sionai  stiffness. 


243  4.  Mode  coupling 

244  We  suggest  that  the  origin  of  the  asymmetrical 

245  shape  of  the  modes  is  caused  by  geometrical 

246  asymmetries,  mostly  due  to  the  tip  positioned  not 

247  exactly  on  the  center  line  of  the  beam.  These 

248  asymmetries  in  the  vibration  pattern  entail  a  cou- 

249  piing  of  flexural  and  torsional  vibrations  via  the  tip 

250  when  the  tip  is  in  force  interaction  with  a  sample 

251  surface.  This  holds  in  particular  when  additional 

252  static  forces  are  applied  in  a  measurement  of 

253  elasticity  [3].  Considering  that  the  tip  with  the 

254  mass  mt  is  offset  from  the  area  center  of  the  cross 

255  section  by  an  amount  d  (see  Fig.  8),  the  flexural 

256  and  torsional  modes  are  no  longer  independent. 

257  They  couple  via  the  forces  acting  on  the  tip  and 

258  hence  lead  to  an  admixture  of  torsional  motion  in 

259  a  flexural  mode  and  vice  versa.  This  coupling  can 

260  be  described  by  the  Eqs.  (1)  and  (8)  taking  account!. 

26 1  the  following  boundary  conditions  at  the  end  with 

262  the  tip  (x  =  L) 

Ely'”  (x,  t)  =  kny(x ,  t)-kj^  Q(x ,  /)  t  ^ 

+  mt (y(x, t)  - ^e(x, t)j  (% 

cTff  (x,  t)  =  —k]h26(x,  t )  -  k„  d Ly(x,  tj  .  k  ' !"- 

+  t n,  d (Ly(x,  t)  -  <j|(x,  /))  (13) 


k\  and  kn  are  the  lateral  and  normal  contact  stiff-  265 
ness,  respectively.  Here,  when  the  cantilever  per-  266 
forms  free  oscillations,  the  boundary  conditions  267 
reduce  to  268 

EI/'(x,  t)  =  m,  (H*, 0  -  j9(x,  ^  (14) 

cT&(x,t)  =mtd(Ly(x,t)  -<||(x(!f||;.  Y  (15) 

This  still  leads  to  a  couphiig^6f?the  modes  because  271 
the  mass  of  the  tip  is  riot 'centered.  Seeking  har-  272 
monic  solutions  y{x^t)  -y(tyeio,t  and  Q(x,  t)  =  273 

0(x)etot  leads  to  the  characteristic  equations  274 


EIy’"(x)  =  mtco2^j^x)  -y{) 


cT  &  (x)  =  mi  d&l(d0(x)  -  Ly{x) )  (17) 

Inserting  Eqs.  (2)  dnd  (10)  into  Eqs.  (16)  and  (17) 
leads  to  the  relation  t]  =  //y2  between  the  flexural 
aridCtopiorial  wave  numbers  y  and  rj  where 


rfis^this  parameter  H  which  determines  the  like¬ 
lihood  of  modal  coupling.  The  mode  shape  of  the 
'ipeak  at  212  kHz  is  a  combination  of  the  flexural 
and  torsional  solution.  The  weighting  of  the 
r  combined  solutions  depends  on  the  type  of  and 
degree  of  coupling  as  well  as  on  the  excitation 
method.  Two  modes  have  frequencies  very  close  to 
each  other.  The  frequencies  of  the  third  flexural 
mode  and  the  first  torsional  mode  lie  within  5%  of 
each  other.  Thus,  during  forced  vibrations,  they 
are  often  excited  simultaneously.  Fig.  9  shows  a 


1 ^ 


AFM  tip 


Fig.  8.  Parameters  entering  in  the  description  of  mode  coupling 
via  the  forces  acting  on  an  AFM  tip  of  length  h  which  is  slightly 
offset  by  an  amount  d  from  the  center-line  of  the  beam. 


Amplitude  distribution  of  the  third  flexural  mode 


Amplitude  distribution  of  the  first  torsional  mode 


Clamped  end 


End  with  the  tip 


Fig.  9.  Spatial  plots  of  modes  3  and  4  calculated  from  the 
flexurai/torsional  coupling  model  with  H  =  0.025.  When  the 
frequencies  are  close  together  small  perturbations  cause  these 
modes  to  couple. 
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292  plot  of  the  combination  of  the  first  torsional  and 

293  the  third  flexural  mode. 

294  A  resonance  was  also  observed  at  265  kHz.  The 

295  corresponding  mode  does  not  fit  into  an  analysis 

296  based  on  mode  coupling  as  explained  above.  Most 

297  likely,  this  mode  originates  from  nonlinear  cou- 

298  pling  into  flexural  motion  about  an  axis  perpen- 

299  dicular  to  the  plan  view.  Such  a  mode  would  have 

300  f\  about  a/b  -  22.5  times  the  usual  fx  for  flexure. 

30 1  This  value  matches  with  the  observed  frequency  at 

302  265  kHz.  Calculating  the  elastic  strain  amplitudes 

303  prevailing  in  the  beams  for  both  modes,  i.e. 

304  Cflcst  ~  {b/2){Sy/6x2)  and  «  0a/£  (Eqs.  (7) 

305  and  (10))  and  using  the  measurements  for  the 

306  amplitudes  as  shown  in  Figs.  (5)  and  (6),  one  ob- 

307  tains  values  of  the  order  £flcx  ^  10“5  and 

308  stors «  10-4.  These  rather  large  values  make  non- 

309  linear  mixing  of  modes  based  on  higher-order 

310  elasticity  possible  [26]. 


311  5.  Conclusion 

312  In  summary,  we  have  shown  that  mode  coupling 

313  due  to  geometrical  and  mass  asymmetries  account 

314  for  a  number  of  resonances  besides  the  standard 

3 1 5  flexural  and  torsional  modes.  The  coupled  modes 

316  can  be  observed  in  many  commercial  cantilevers. 

317  They  are  sometimes  called  spurious  modes,  yet 

318  they  can  be  explained  using  quite  simple  physical 

319  principles.  The  large  strains  values  obtained  ih 

320  dynamic  operations  may  lead  to  nonlinear  mixing 

321  of  modes. 
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Abstract 

Atomic  force  acoustic  microscopy  (AFAM)  can  be  used  to 
make  quantitative  images  of  surface  stiffness  with  high  resolu¬ 
tion.  The  technique  utilizes  the  dynamic  response  of  the  atomic 
force  microscope  (AFM)  cantilever,  specifically  in  terms  of  the 
higher-order  cantilever  modes.  Motion  of  the  specimen  while  in 
contact  with  the  AFM  tip  induces  vibrations  in  the  cantilever. 
The  dynamical  properties  of  the  cantilever  change  during  scan¬ 
ning  due  to  the  variations  of  surface  stiffness.  Both  stiff  and  com¬ 
pliant  cantilevers  may  be  used  for  stiffness  imaging  if  the  appro¬ 
priate  vibration  mode  is  chosen.  Experimental  AFAM  results  on 
polysilicon  for  microelectromechanical  systems  (MEMS)  are  pre¬ 
sented.  The  AFAM  results  provide  stiffness  information  about 
the  material  microstructure  that  is  not  available  from  topogra¬ 
phy  information  alone.  The  technique  provides  in-situ  stiffness 
information  without  damaging  the  MEMS. 

Introduction 

The  atomic  force  microscope  (AFM)  was  developed  as  an  in¬ 
strument  for  measuring  surface  topography  while  operating  in  a 
quasistatic  mode  [1].  The  deflection  of  the  microscope  cantilever 
plotted  as  a  function  of  surface  location  gives  a  high-resolution 
image  of  surface  topography.  More  recent  applications  of  the 
AFM  have  exploited  the  dynamic  characteristics  of  the  AFM. 
Many  versions  of  dynamic  atomic  force  microscopy  have  been 
proposed  for  imaging  specimens  and  measuring  material  proper¬ 
ties  [2-5].  All  of  these  methods  rely  on  the  relative  motion  be¬ 
tween  the  atomic  force  microscope  (AFM)  tip  and  the  specimen 
.  surface.  The  high  frequency  excitation  of  the  specimen  surface 
or  AFM  cantilever  can  provide  a  better  signal-to-noise  ratio  for 
imaging  while  doing  less  damage  to  the  specimen  surface.  These 
techniques  utilize  the  dynamic  response  of  the  cantilever,  specifi¬ 
cally  in  terms  of  the  vibrational  modes.  The  dynamic  response  of 


an  elastic  cantilevered  beam  in  contact  with  a  surface,  however, 
is  not  easily  described.  The  tip-sample  interaction  forces  are' 
in  general,  nonlinear  and  'fixe  the  subject  of  intensive  research 
themselves.  When  an  elastic  beam  interacts  with  these  surface 
forces,  the  many  flexural  and  torsional  modes  must  be  included 
as  part  of  the  complete  dynamics.  These  higher  modes  are  often 
excited  simultaneously  in  experiments. 

In  much  of  the  current  research,  the  amplitude  of  the  sur¬ 
face  motion  is  kept  small.  In  this  case,  the  motion  may  be 
studied  by  linearizing  the  nonlinear  tip-sample  interaction  forces 
about  an  equilibrium  point.  The  cantilever  motion  may  then 
be  modeled  as  a  beam  with  a  linear  spring  attached  near  the 
end.  Measured  natural  frequencies  of  the  cantilever  with  at¬ 
tached  spring  then  define  the  spring  stiffness.  Appropriate  tip- 
sample  interaction  models  are  then  used  to  relate  the  spring 
stiffness  to  the  Young’s  modulus  of  polysilicon.  Atomic  force 
acoustic  microscopy  (AFAM)  is  anticipated  to  have  many  uses 
for  in  situ  measurements  of  polysilicon  microelectromechanical 
systems  (MEMS). 


Atomic  Force  Acoustic  Microscopy 

The  atomic  force  microscope  cantilever,  as  shown  schemati¬ 
cally  in  Fig.  1,  is  modeled  as  an  Euler-Bemoulli  beam.  The 
beam  is  assumed  uniform  and  homogenous  with  constant  cross- 
section.  It  is  cantilevered  at  one  end.  At  the  end  opposite  of 
the  cantilever  a  tip  with  small  radius  is  attached.  The  sample  is 
offset  such  that  the  beam  has  an  initial  static  deformation  y  (a;) . 
Motion  about  this  initial  equilibrium  deformation  is  given  by 
w  (x,  t) .  The  linearized  boundary  value  problem  for  this  system 
has  a  governing  equation  and  boundary  conditions  given  by  [6,7] 


dx*  +P  dft  ~U’ 


(1) 


79 


w  (x  =  0,  t)  =  0, 


vJ  (x  =  0,t)  =  0, 


(2) 


Tip  Contact 


w"  (x  —  L, 1)  =  0,  (3) 

JETu/"  (s  =  L,t)  =  /«/;(*  =  L,t).  (4) 

In  Eqs.  (l)-(4),  w(x,  t)  defines  the  dynamic  cantilever  position 
relative  to  y  (a:).  The  AFM  cantilever  is  defined  by  the  modulus 
E}  the  area  moment  of  inertia  I,  the  volume  density  p,  and  the 
cross-sectional  area  A,  Here,  El  and  pA  are  assumed  uniform 
over  the  length  of  the  cantilever.  The  boundary  conditions  given 
by  Eqs.  (2)  correspond  to  conditions  of  zero  displacement  and 
zero  slope  at  x  =  0.  The  boundary  condition  given  by  Eq.  (3) 
corresponds  to  zero  moment  at  x  =  L. 

The  force  balance  between  the  shear  force  of  the  beam  and 
the  force  from  the  contact  mechanics  is  given  by  Eq.  (4).  The 
contact  force  is,  in  general,  nonlinear.  Here,  the  motion  is  as¬ 
sumed  small  such  that  a  linear  contact  force  may  be  used.  The 
contact  force  is  proportional  to  the  displacement  of  the  end  of 
the  beam  with  the  effective  contact  stiffness  k  (linear  spring  con¬ 
stant)  as  the  proportionality  constant.  Because  a  linear  model  is 
used  to  describe  the  tip-sample  interaction  force,  the  results  are 
restricted  to  small  tip  displacements.  Such  a  condition  is  easily 
satisfied  experimentally.  The  actual  value  of  k  is  given  in  terms 
of  the  material  properties  of  the  tip  and  sample  and  depends  on 
the  interaction  model  used.  Simple  Hertzian  theory  [8-9],  more 
complex  Maugis  mechanics  [9],  or  other  models  [10-12]  may  all  be 
linearized  about  some  initial  equilibrium  position  to  determine 
this  spring  constant. 

Eqs.  (l)-(4)  define  completely  the  linearized  flexural  vibration 
problem.  The  natural  frequencies  of  the  cantilever  vibrations  are 
dependent  on  the  linear  spring  constant  k.  Solution  of  the  corre¬ 
sponding  eigenvalue  problem  leads  to  the  characteristic  equation 
that  governs  the  allowable  frequencies.  The  characteristic  equa¬ 
tion  is  given  by  [6] 


p 


w&t) 


AFM  Beam 


I  Offset,  zq 
ISample  | 

‘  Static  deflection,  y(x) 


Figure  1:  Schematic  of  the  problem.  The  AFM  probe  is  in  con¬ 
tact  with  a  specimen.  Initial  contact  is  made  when  the  sample 
offset,  zq,  is  zero.  The  static  sample  offset  causes  a  static  beam 
deflection,  y(x).  The  dynamic  motion,  w(x,t),  is  defined  relative 
to  y(x).  The  contact  forces  between  the  tip  and  sample  are  may 
be  discussed  for  different  contact  models. 


Contact  Mechanics 

The  mechanics  of  the  AFM  tip  as  it  interacts  with  a  speci¬ 
men  are  complex  and  generally  nonlinear.  Here,  the  contact  is 
assumed  frictionless  and  the  contact  force  is  assumed  to  act  nor¬ 
mal  to  the  tip.  The  geometry  of  the  contact  problem  is  shown  in 
Fig.  (2).  The  goal  is  to  relate  the  applied  force  p  to  the  penetra¬ 
tion  depth  of  the  tip  into  the  specimen  S.  Important  parameters 
for  this  problem  are  the  radius  of  the  contact  area  a,  the  work 
of  adhesion  at  contact  u,  and  the  stress  in  the  cohesive  zone  <r0. 

The  material  properties  and  geometry  of  the  specimen  and 
AFM  tip  also  are  important.  These  parameters  are  shown  in 
Fig.  (3).  The  reduced  tip-sample  curvature, 


R 


Rt 


(7) 


accounts  for  curvature  of  both  tip  (Rt)  and  specimen  (Ra).  Most 
often,  the  specimen  is  assumed  flat  such  that  R  =  Rt .  The  re¬ 
duced  elastic  modulus  K  is  defined  by 


73  (cos7cosh7  + 1) 

-  gjTjJp  (sinh7Cos7  -  sin7Cosh7)  =  0.  (5) 

Equation  (5)  is  solved  numerically  for  the  values  of  7.  The  dis¬ 
persion  relation 

is  then  used  to  determine  the  natural  frequencies  in  terms  of  the 
wave  numbers.  The  linear  vibrations  of  the  AFM  may  be  used 
to  study  the  material  properties  of  specimens.  Knowledge  of  the 
natural  frequencies  can  be  used  to  determine  the  effective  spring 
constant  k.  The  tip-sample  mechanics  are  then  used  to  deduce 
the  local  Young’s  modulus.  The  most  general  contact  mechanics 
model,  due  to  Maugis  [9]  is  discussed  below.  The  Maugis  model 
was  first  introduced  to  the  AFM  literature  by  Burnham  and 
coworkers  [13]. 


1  ,  i-*\ 

K  4  \Et  Ea  )’ 


(8) 


where  Et ,  vtf  E9)  and  us  are  the  Young’s  modulus  and  Poisson’s 
ratio  of  the  rip  and  sample,  respectively.  The  general  goal  for 
materials  characterization  studies  with  AFAM  is  to  extract  and 
image  Young’s  modulus  of  the  specimen.  Such  a  procedure  in¬ 
volves  measurement  of  K.  The  material  properties  of  the  tip  are 
generally  assumed  to  be  known.  In  addition,  a  calibration  proce¬ 
dure  or  scanning  electron  microscopy  is  necessary  to  determine 
the  tip  radius. 

The  most  comprehensive  tip  mechanics  model  is  that  of 
Maugis  [9,13,14].  The  Maugis  model  captures  many  aspects  of 
the  contact  mechanics  observed  in  AFM  studies  through  the  use 
of  fracture  mechanics  concepts.  Most  importantly,  it  can  account 
for  adhesion  hysteresis.  The  Maugis  model  provides  a  continuous 
functional  dependence  from  Hertzian  (or  DMT  [11])  mechanics, 
without  adhesion  hysteresis,  to  the  JKR  [12]  limit  of  maximum 
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Figure  2:  Schematic  of  the  parameters  of  the  contact  mechanics. 
An  applied  load  p  results  in  a  penetration  to  depth  6 .  The  contact 
radius  is  given  by  a  and  the  radius  to  the  edge  of  the  cohesive 
zone  is  denoted  by  c. 


adhesion  hysteresis.  The  dimensionless  parameters  used  are  de¬ 
fined  as 


C 

I 


f 1 


( 

1 

f 

( 

i 


M=*)-  « 

A  -  ‘{SsT*’  (>») 

where  A  is  the  contact  radius,  P  is  the  applied  force,  A  is  the 
penetration  depth,  and  A  is  the  adhesion  parameter. 

In  terms  of  these  parameters,  several  different  contact  me¬ 
chanics  models  may  be  described.  Several  of  these  models  are 
illustrated  by  the  force-penetration  curves  (P  vs.  A)  shown  in 
Fig.  (4).  Hertzian  mechanics  is  given  by 


pH  =  a3,  Ah  =  A2,  (11) 

which  may  be  written  PH  =  A^2.  In  the  Hertz  model,  adhesion 
is  neglected  such  that  there  is  no  attractive  force  associated  with 
loss  of  contact.  The  DMT  model  is  a  modified  Hertzian  model 
[11].  In  this  case,  a  constant  attractive  force  is  added  to  the 
Hertzian  mechanics,  p£>  =  A3  —  2.  However,  the  tip  geometry  is 
assumed  not  to  change,  AD  =  A2.  The  DMT  PD  (A)  is  noted 
in  Fig.  4  for  A  =  0.  The  Hertz  and  DMT  models  do  not  include 
adhesion  hysteresis  effects.  Thus,  they  are  most  applicable  to 
AFM  research  for  low  adhesion  specimens,  hard  materials  or 
small  radii  of  curvature. 


Figure  3:  Geometry  of  the  contact  problem  in  terms  of  the  tip 
parameters  (subscript  t)  and  specimen  parameters  (subscript  s) . 


In  the  high  adhesion  limit,  the  JKR  model  is  appropriate.  It 
is  given  by  [12] 

Pj  =  A3  -  AV6A:  Aj  =  A2  —  |V6A.  (12) 

o 

This  model  is  most  applicable  for  high  adhesion  specimens,  soft 
materials  or  large  radii  of  curvature. 

The  model  developed  by  Maugis  is  comprehensive  in  that  it 
captures  all  types  of  tip-sample  mechanics  through  the  introduc¬ 
tion  of  an  adhesion  parameter  A.  Attractive  (tensile)  stresses  are 
assumed  to  act  between  the  tip  and  specimen  in  a  cohesive  zone 
that  surrounds  the  contact  area.  These  stresses  are  assumed  con¬ 
stant  (<tq)  within  an  annulus  of  width  d.  For  the  Maugis  model, 
the  force  and  deflection  are 


Pm  =  A3  -  XA2Mu  (13) 

A  m  =  A2  —  ^AX\/ m2  —  1, 


where  _ 

M1  =  /m2  -  1  +  m2  tan"1  \/m2  - 1.  (14) 


In  Eqs.  (13)  and  (14),  m  —  c/a  is  the  ratio  of  the  radius  to 
the  edge  of  the  cohesive  zone  to  the  radius  of  the  contact  area. 
The  width  of  the  cohesive  zone  d  =  c—a.  An  additional  equation 
is  necessary  for  the  Maugis  mechanics.  The  Griffith  equilibrium 
equation  defines  the  relation  between  A,  m,  and  A.  It  is  given  bv 

[9] 


AA2  4A2Aw 
——M2  H - —  M3  -  1, 


(15) 
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where 


M2  =  \/m2  —  1  +  (m2  —  2)  tan  1  y/rrfi  —  1,  (16) 

M3  =  1  —  \/m2  -  1  tan"1  y/m1—  1.  (17) 

A  comparison  of  this  contact  model  for  several  different  val¬ 
ues  of  A  is  shown  in  Fig.  4.  When  A  is  positive,  the  sample 
is  deflected  away  from  the  AFM  tip.  A  negative  value  for  A 
implies  that  the  material  is  bulging  out  from  the  surface.  The 
Maugis  model  includes  attractive  forces  from  adhesion  such  that 
the  cantilever  is  deflected  downward  when  contact  is  initiated.  In 
this  case,  the  tensile  stresses  in  the  cohesive  annulus  are  greater 
than  the  compressive  stresses  in  the  contact  area.  The  effect  of 
adhesion  hysteresis  is  also  observed  in  Fig.  4.  Several  curves 
shown  in  Fig.  4  bend  back  on  themselves.  These  curves  have 
a  section  of  positive  slope  when  both  P  and  A  are  negative. 
Points  on  this  part  of  the  curve  are  unstable  equilibria.  Thus, 
when  moving  continuously  on  these  curves  loss  of  contact  occurs 
when  the  cantilever  stiffness  exceeds  the  contact  stiffness.  In  ad¬ 
dition,  when  contact  is  reinitiated,  the  tensile  cohesive  stresses 
would  pull  the  AFM  tip  back  into  contact.  This  type  of  hys- 
teretic  behavior  results  for  all  values  of  A  that  are  larger  than 
about  0.94.  The  Maugis  model  reduces  to  that  of  DMT  mechan¬ 
ics  when  A  =  0.  In  the  opposite  high  adhesion  limit,  when  A  is 
large,  the  JKR  limit  [12]  is  approached. 


Figure  4:  Normalized  force-displacement  relation  for  various  val¬ 
ues  of  A  using  the  Maugis  contact  mechanics  model. 

Surface  Stiffness 

As  the  AFM  tip  is  placed  in  contact  with  the  specimen,  it  ini¬ 
tially  settles  at  an  equilibrium  position  somewhere  on  the  P  (A) 
curve.  The  AFM  cantilever  will  have  some  initial  deflection  y  (#) 


determined  by  the  force  balance  at  the  tip.  The  cantilever  or 
sample  is  then  vibrated  at  small  amplitude.  If  the  displacements 
of  the  AFM  cantilever  remain  small,  a  linearization  of  the  above 
contact  models  may  be  used  to  determine  the  appropriate  spring 
stiffness  necessary  for  the  vibration  problem.  In  this  case,  the 
linear  spring  constant  (in  nondimensional  form)  is  defined  as 


In  dimensional  form,  we  have 


1/3  dP 
dA' 


For  the  Maugis  model,  the  relation  of  the  surface  stiffness  is 
found  from  Eqs.  (13)-(15).  The  force  and  deflection  are  both 
given  in  terms  of  the  contact  area  and  parameter  A.  In  addition, 
the  Griffith  equilibrium  equation  shows  that  A  is  also  a  function 
of  A.  Thus,  we  have 


=  3A2-M1  ^2XA  +  A2^  , 


•  (21) 


The  rate  of  change  of  A  with  A  is  found  from  Eq.  (15)  to  be 

d\  2 A  3AM2+'4AAf3  ,  . 

dA  A  3AM%  H-  I6AM3  ’  {  J 

where  M2  and  M3  are  defined  in  Eqs.  (16)  and  (17).  Thus,  the 
contact  stiffness  from  the  Maugis  model  is 

€u  =  5^2,  (23) 


Q  =  3 A2M2  +  I6AXM3  -  8X2M1M3, 


R  =  ZM2A  (3A  +  2Av/(m2  -  1)) 

+  I6M3A  (3A  -  Ai/(m2  -  1))  .  (25) 


The  general  Maugis  result  for  the  stiffness  reduces  to  the  Hertz 
and  DMT  cases  (both  have  the  same  £  vs.  P  curves)  for  A  — *  0 
and  to  the  JKR  result  when  A  is  very  large.  Figure  (5)  is  a  plot 
of  £  as  a  function  of  applied  load  P  for  various  values  of  A. 


Experiments 

The  experiment  used  here  is  similar  to  that  described  else¬ 
where  [2].  A  standard  atomic  force  microscope  (Thermomicro¬ 
scopes,  Autoprobe  CP)  is  modified  for  AFAM  capabilities.  A 
schematic  of  the  AFAM  experiment  is  shown  in  Fig.  (6).  The 
specimen  is  placed  on  a  broadband  ultrasonic  transducer  (2.25 
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Figure  5:  Normalized  surface  stiffness  as  a  function  of  applied 
load  P,  for  various  values  of  A  using  Maugis  contact  mechanics. 

MHz  center  frequency).  The  transducer  is  mounted  on  the  AFM 
stage  and  is  excited  by  the  function  generator.  The  output  from 
the  AFM  photodiode  is  mixed  in  the  lock-in  amplifier  with  the 
function  generator  output.  The  lock-in  output  and  the  photodi¬ 
ode  output  are  both  used  by  the  AFM  computer  for  imaging. 


Figure  6:  Schematic  of  the  atomic  force  acoustic  microscopy  ex¬ 
periment.  Both  the  topography  and  lock-in  amplifier  output  are 
used  for  imaging. 

Two  sets  of  AFAM  results  axe  presented  in  Figs.  (7)-(10). 
The  specimens  are  polysilicon  strips  fabricated  at  the  CRONOS 
MEMS  foundry.  The  scan  size  of  each  figure  is  2  //m  x  2  fzm. 

Topography  and  AFAM  results  using  mode  1  of  a  stiff  can¬ 
tilever  (40  N/m)  are  shown  in  Figs.  (7)  and  (8),  respectively. 
The  applied  force  was  1000  nN.  The  lightest  color  in  the  topog¬ 
raphy  image  is  85  nm  higher  than  the  darkest  color.  The  topog¬ 
raphy  image  clearly  shows  a  few  of  the  grains  of  the  polysilicon. 
The  AFAM  image,  shown  in  Fig.  (8),  was  made  using  an  excita¬ 


tion  frequency  of  816  kHz.  The  contrast  in  this  image  is  related 
to  shifts  in  the  cantilever  resonance  peak.  The  chosen  excita¬ 
tion  frequency  was  slightly  higher  than  the  resonance  frequency. 
Thus,  lighter  colors  implies  less  stiffness.  The  darkest  colors  are 
difficult  to  interpret.  The  peak  may  have  shifted  greatly  higher 
or  lower.  The  stiffness  contrast  (AFAM  image)  of  the  polysilicon 
is  considerable.  Within  individual  grains,  contrast  variations  are 
observed.  In  particular,  the  grain  indicated  in  Fig.  (7)  with  the 
arrow  is  especially  interesting.  The  AFAM  image  indicates  that 
this  grain  has  large  variations  in  stiffness  throughout.  The  grains 
are  expected  to  be  single  crystal  grains  which  grow  perpendicu¬ 
lar  to  the  MEMS  surface.  The  source  of  these  stiffness  variations 
is  still  under  investigation. 


Figure  7:  Topography  image  of  polysilicon  using  a  stiff  cantilever 
(40  N/m).  The  scan  size  is  2  /xm  x  2  fim 

Topography  and  AFAM  results  using  a  weaker  cantilever  (3 
N/m)  are  shown  in  Figs.  (9)  and  (10),  respectively  at  a  different 
location  on  the  specimen.  The  topography  image,  Fig.  (9),  has 
a  maximum  height  difference  of  85  nm.  The  lightest  color  corre¬ 
sponds  to  the  highest  topography.  The  AFAM  image  is  shown  in 
Fig.  (10).  These  images  were  made  using  the  second  vibration 
mode  of  a  compliant  cantilever  with  natural  frequency  in  contact 
of  623  kHz.  In  this  image,  the  lightest  color  corresponds  to  an 
amplitude  of  1.2  mV.  The  applied  force  on  the  cantilever  was  700 
nN.  Again  the  stiffness  variations  within  the  grains  are  clearly 
observed.  .  - 

Effective  Cantilever  Length 

AFM  cantilevers  are  normally  manufactured  using  etching 
processes.  These  techniques  are  not  currently  precise  enough 
for  the  cantilever  geometry  to  be  known  with  high  accuracy.  For 
AFAM  measurements  the  placement  of  the  tip  along  the  can¬ 
tilever  is  a  critical  piece  of  information.  SEM  images  show  that 


Figure  8:  AFAM  image  using  the  first  cantilever  mode  of  a  stiff 
cantilever.  The  first  free  natural  frequency  is  170  kHz.  The  scan 
size  is  2  pm  x  2  pm. 

the  tip  may  be  anywhere  from  9099%  of  the  cantilever  length 
from  the  fixed  end.  Attempts  to  characterize  the  effective  length 
of  the  cantilever  (i.e.,  the  distance  from  the  fixed  end  to  the  tip 
location)  have  met  with  some  success  [2].  Here  an  alternative 
method  is  proposed. 

The  shift  in  the  natural  frequencies  of  the  beam  are  affected 
by  the  contact  stiffness.  This  stiffness  is  a  function  of  the  applied 
load.  If  the  Hertzian  contact  model  is  used  we  have 

k=^K^RP.  (26) 

For  a  certain  applied  load  p\  and  values  of  R  and  K  there  is 
a  unique  corresponding  stiffness  Ky.  Now  imagine  that  R  and  K 
are  held  constant  while  another  load  P2  is  applied  with  corre¬ 
sponding  stiffness  «2-  If  P2  =  Bpu  where  B  is  a  constant  factor, 
then  K2/K1  =  B 1^3.  Thus,  the  shift  in  stiffness  as  a  function  of 
relative  applied  load  is  known.  Also,  the  shifted  frequencies  cor¬ 
responding  to  the  different  values  of  k  may  be  determined.  For  a 
set  of  applied  forces,  there  is  one  effective  length  which  provides 
the  best  fit  to  the  data  for  a  single  value  of  L.  To  illustrate  the 
method  consider  the  following.  A  cantilever  with  first  free  fre¬ 
quency  of  176  kHz  is  placed  in  contact  with  Si  <100>  with  an 
applied  load  of  800  nN.  The  first  natural  frequency  is  measured 
to  be  844  kHz  as  shown  In  Fig.  (11).  This  frequency  is  4.8  times 
larger  than  (/1)  jree .  If  the  cantilever  were  pinned  at  the  end,  the 
maximum  shift  possible  is  3.92  times  greater  than  the  free  [15]. 
Thus,  the  effective  length  has  shortened  such  that  the  frequency 
is  higher.  Next,  another  higher  load  is  applied,  such  as  1600  nN. 
The  contact  stiffness  should  increase  by  a  factor  of  21/3  with  a 
known  change  in  frequency.  Other  loads  are  also  applied  in  a 


Figure  9:  Topography  image  of  polysilicon  using  a  compliant 
cantilever  (3  N/m).  The  scan  size  is  2  pm  x  2  pm 

similar  manner.  Finally,  .the  experimental  results  are  compared 
with  curves  of  natural  frequencies  for  different  values  of  the  ef¬ 
fective  length.  The  curve  for  Leff  =  0.969X  was  found  to  be  the 
best  fit  to  these  data.  This  result  is  the  unique  solution  for  both 
the  effective  length  and  contact  stiffness.  It  should  be  noted  that 
this  method  requires  no  information  about  the  tip  radius  R  nor 
the  material  properties  of  the  beam.  The  required  measurements 
are  relative  to  the  free  experimental  measurements.  Such  a  tech¬ 
nique  may  be  easily  extended  to  non-uniform  cross-sections.  The 
effect  of  damping  is  expected  to  have  an  influence  on  this  result, 
but  was  not  considered.  Comparison  of  this  method  for  deter¬ 
mining  the  effective  length  with  other  methods,  such  as  SEM, 
are  currently  in  progress.  The  technique  is  restricted  to  stiff  can¬ 
tilevers  with  high  applied  loads  such  that  the  Hertzian  theory  is 
valid. 

Young’s  Modulus  for  Pofysilicon 

Estimates  of  Young’s  modulus  for  polysilicon  are  now  made 
based  on  the  results  presented  above.  Both  the  first  mode  of 
the  stiff  cantilever  (40  N/m)  and  the  second  mode  of  the  more 
compliant  cantilever  (3  N/m)  are  used.  For  the  stiff  cantilever 
the  effective  length,  Leff  =  0.969L  was  used  to  determine  the 
contact  stiffness  on  the  polysilicon  MEMS  material.  Again  the 
relative  frequency  shift  is  examined.  If  the  applied  load  is  the 
same  on  both  the  silicon  and  on  the  poly  silicon,  the  ratio  of  the 
contact  stiffnesses  is 


Equation  (27)  is  rearranged  and  solved  for  Epoiy  where  use  of 
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Figure  10:  AFAM  image  using  the  second  cantilever  mode  of  a 
compliant  cantilever  (623  kHz).  The  first  free  natural  frequency 
is  90  kHz.  The  scan  size  is  2  /im  x  2  pm. 


the  reduced  modulus,  Eq.  (8)  has  been  used.  The  result  is 


E+pcly  : 


[1  -  F]  ^ 


(28) 


£si 


where  the  factor  F  is  given  by 


The  reduced  radius  on  silicon  Bs\  is  assumed  to  be  that  of  the 
tip  radius  Rt  because  the  silicon  specimens  were  very  fiat.  The 
reduced  radius  for  the  polysilicon  must  reflect  the  topography 
seen  in  Figs.  (7)  and  (9).  In  this  case, 


{r&v  ^  (30) 

where  Rs  is  the  sample  radius.  The  sample  radius  is  estimated 
to  be  on  the  order  of  50  nm.  The  tip  radius,  measured  using 
calibration  gratings,  was  Rt  =  120  nm.  Using  vt  =  0.181,  Et  = 
130.4  GPa,  and  i/^iy  =  0.18,  gives  the  value  for  Epoiv  «  104 
GPa .  This  value  for  Ep^y  Is  about  65%  of  that  found  from  tensile 
tests  of  polysilicon  [16].  The  present  results  are  encouraging,  but 
still  need  improvement. 

For  the  more  compliant  cantilever,  the  adhesion  effects  are 
much  more  important.  Therefore,  the  method  used  for  the  stiff 
cantilever  is  not  expected  to  apply.  The  tip  radius  was  mea¬ 
sured  using  calibrated  gratings  and  determined  to  be  Rt  =  160 
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Figure  11:  Experimental  frequecy  spectra  for  three  different  ap¬ 
plied  loads.  The  free  natural  frequency  of  the  cantilever  (40 
N/m)  was  176  kHz.  The  specimen  was  silicon  <100>. 


nm.  The  applied  load  was  700  nN.  The  measured  second  natural 
frequency  of  623  kHz  corresponds  to  a  contact  stiffness  of  k  — 
41.1  N/m.  The  measured  pull-off  force  was  112  nN.  From  the 
-definition  of  the  dimensionless  force,  Eq.  (13),  this  pull-off  force 
implies  that  P  «  6.25  for  the  applied  load  used.  From  Eq.  (23), 
the  range  of  values  for  £  =  dP/dA  «  2.9  -  3.0  for  all  values  of  A. 
Rearranging  Eq.  (19)  gives 


K  = 


K 


3/2 


\{^R)R\^^ 


(31) 


Substitution  of  the  above  values  into  Eq.  (31)  gives  K  = 
400  MPa .  In  terms  of  the  specimen  stiffness  we  have 


Ea  = 


1-2/2 

4  JL  ,  1  -Vi  4 
3  Bt 


(32) 


Using  ut‘=  0.181,  Et  =  130.4  GPa,  and  us  =  0.18  gives 
Es  —  290 MPa.  This  value  is  unrealistically  low  for  polysili¬ 
con.  This  result  highlights  the  difficulties  associated  with  mak¬ 
ing  measurements  of  stiff  materials  with  compliant  cantilevers. 
Sources  of  error  associated  with  these  measurements  include  the 
tip  radius,  cantilever  stiffness,  and  pull-off  force.  Errors  associ¬ 
ated  with  pull-off  tests  for  measuring  adhesion  have  been  noted 
by  others  [14]. 
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Conclusions 

In  this  article,  the  use  of  atomic  force  acoustic  microscopy 
(AFAM)  has  been  discussed  as  applied  to  polysilicon  micro- 
electromechanical  systems  (MEMS).  An  estimate  of  Young’s 
modulus  of  polysilicon  has  been  made  with  this  technique.  The 
results  presented  here  are,  unfortunately,  much  lower  than  pre¬ 
vious  direct  measurements.  However,  several  improvements  of 
the  measurement  technique  may  prove  useful.  The  dependence 
on  tip  radius  is  substantial.  The  calibrated  grating  method  for 
radius  measurement  has  not  yet  been  validated  as  an  effective 
method  to  determine  R.  In  addition,  effects  from  surface  forces 
may  substantially  alter  the  effective  value  for  the  contact  load 
P.  The  value  used  in  the  above  calculations  corresponds  to  all 
forces  sensed  by  the  cantilever.  All  of  these  sources  of  error  are 
being  examined  in  order  to  improve  the  measurements.  The  use 
of  AFAM  for  other  materials  has  shown  the  potential  of  this 
technique  for  in  situ  measurements  of  MEMS. 
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added  to  the  numerical  model.  Viscous  damping  was  included  in  the  form  of  a  dashpot  with 


damping  constant  c  in  parallel  with  the  spring  corresponding  to  k*.  For  the  minimization  in  Eq.  (1).  Each  entry  in  Table  II  represents  the  average  of  four  data  sets;  as  described  in 

procedure,  Li/L  and  k*  were  first  varied  to  minimize  the  error.  Then,  an  appropriate  value  Sect.  II  A,  each  data  set  typically  yielded  six  values  of  Mm.  The  uncertainties  quoted  in 

of  c  was  added  to  further  reduce  the  error.  The  process  was  repeated  until  the  required  error  Table  II  represent  the  standard  deviation  of  these  multiple  measurements.  The  accuracy  of 

level  was  achieved.  For  this  cantilever,  the  range  over  which  L\jL  was  allowed  to  vary  was  the  value  used  for  Mrej  is  potentially  an  additional  source  of  measurement  uncertainty.  If  the 
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TABLE  IIL  Values  for  the  indentation  modulus  M  of  Nb  (in  GPa).  A  range  of  values  for  bulk 
Nb  obtained  from  the  literature  is  shown.  Also  included  are  results  obtained  by  surface  acoustic 
wave  spectroscopy  (SAWS)  and  instrumented  indentation  techniques  (IIT)  on  the  same  sample 
measured  by  AFAM.  The  minimum  and  maximum  AFAM  values  from  Table  II  are  eiven. 


The  nonlinear  analysis  which  follows  requires  a  more  amenable  form  for  the  nonlinear  contact. 


with  boundary  conditions  problem  as 


subscripts  on  An  and  are  implicit  throughout  the  remainder  of  the  article.  The  linear  mode 


where 


Next  we  define  ip  —  aTi  —  q  such  that  qf  =  ip'  —  a.  Separating  Eq,  (56)  into  real  and  imaginary 
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5.1  Characteristic  Equation. 


combinations  of  AFM  and  ultrasonics  have  been  proposed  [9] [14] [15]. 

The  most  significant  difference  of  the  dynamic  methods  versus  the  static 


Figure  2.2.  Ultralever  chips  with  four  cantilevers  [28,  pp.  4-32]. 


It  can  be  seen  from  Figure  (2.7)  that  Hertz,  DMT,  BCP  are  not  hysteretic,  whereas 
JKRSis. 


sample,  the  surface  forces  will  become  negligible.  Under  such  circumstances  the  Hertzian 


contact  mechanics  theories,  the  Hertzian  theory  is  the  simplest  with  the  assumption 


to  non-conservative  forces. 


where  y-  kL  is  the  normalized  wave  number.  From  this  equation  the  allowable  wave  observed  from  this  figure,  the  resonance  frequencies  acquired  from  the  characteristic 

numbers  can  be  found.  U.  Rabe,  K.  Janser  and  W.  Arnold  have  done  extensive  work  on  equation  (3.5.2)  show  a  consistent  trend  of  increasing  with  ratio  {JOkc).  When  K/kc 


contact  boundary  condition.  However,  they  used  the  Johnson-Kendall-Roberts  model 


Note  that  these  spring  constants  are  not  independent.  They  are  all  defined  by  the  material 
properties  of  the  tip  and  specimen  together  with  the  initial  offset.  The  nonlinear  constants 


Substituting  equations  (4.2.11)  and  (4.2.12)  into  equation  (4.2. 13)  gives 


encompasses  a  continuous  variation  of  the  adhesion  parameter  X.  In  the  limit  of  m  going 


to  1  +  £,  the  Maugis  theory  will  approach  the  JKRS  theory.  Note  that 


of  the  root  Of  (4.3.16).  Taylor’s  scries  expansion  of  the  contact  boundary  condition  holds  true  only  if  w(L,t)  <  A 

as  shown  in  (3.4.4).  Thus  a\  must  satisfy  a\  <  A. 

Since  the  form  of  the  force-displacement  curve  is  concave,  the  beam  vibrations  are  The  vibration  amplitude  a\  and  the  resonance  frequencies  of  the  beam  are  closely 

expected  to  act  according  to  nonlinear  softening.  Thus,  from  the  physical  perspective,  the  related.  The  characteristic  equation  is  solved  numerically  and  the  amplitude-frequency 


safely  applied  to  other  nonlinear  systems  with  z  less  than  0. 1. 


ascending  order,  is  fin,  fix,  fi\A,  fix,  fix,  fin,  fix,  fin,  fiiS,  fii5.  After  the  intersection  point 
fiAs,  the  sequence  of  sensitivity  level  for  the  first  five  modes  is  completely  reversed  and 


difficulties  in  experiments  for  accurate  measurements  of  the  frequency  shift. 


In  the  Figures  (5. 1)  and  (5.2),  J5  is  100,  and  p  is  varied  from  1  to  100.  For  most  AFM  than  the  other  modes,  especially  at  p  near  7,  where  the  damping  of  the  second  mode 

measurements,  p  is  less  than  100.  It  can  be  seen  that  the  first  and  second  mode  curves  reaches  maximum. 


increases  quickly.  The  imaginary  parts  of  the  other  inodes  are  well  below  that  of  the 


imaginary  part  of  the  second  mode  with  respect  to  p  and  small  variation  in  /?. 


inodes  intersect  are  different:  in  Figure  (5. 1),  it  is  about  7.3;  in  Figure  (5.6),  it  is  about  8; 


more  than  those  of  others  when  amplitude  increases.  In  Figure  (5.9),  the  influence  of 


variables  are  introduced  according  to  Tn  =  ent,  n  =  0,12 . each  of  which  represents  a 


Similarly,  by  combining  all  the  coefficients  of  the  £*  order  terms,  the  £*  order  equation 


with  the  same  boundary  conditions  as  (6.2.4). 


From  equation  (6.4.10),  fyois  substituted  into  equations  (6.4.11)  to  yield 


r  ■  ,  co2 -Cl2  (o)2- Cl2 )2  J  K  h  system  are  easier  to  measure. 

The  real  part  of  7}r(t)  is  For  m  undamPed  linear  system,  if  the  frequency  of  the  driving  force  is  very  close  to 

the  primary  natural  frequency  ( Ob)  of  the  system,  the  oscillation  of  the  system  will 


become  unbounded  regardless  of  the  amplitude  of  the  excitation.  In  a  continuum  system 


(»)  3pnji]<fcuv 


(»)  opmFlduiv 


damping  on  the  response  curve. 


amplitude  is  inversely  proportional  to  damping.  Figure  (6.7)  shows  the  influence  of 


Figure  6.10.  Response  curve  of  the  fourth  mode  for  /?=50,  c  -  0.01. 


frequency.  By  the  method  of  harmonic  balance,  the  solution  is  written  as 


r}r^Ar(Tx)eJ^^hrem\  A,=  — .  (6.4.51) 


Then,  by  the  characteristic  equation. 
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Scientific  Instruments ,  Vol.  71,  No.  6, 2403-2408  (2000) 


The  most  important  component  of  the  AFM  is  the  cantilever  used  to  probe  the 
sample  surface.  The  cantilevers  tend  to  be  very  small  and  are  fabricated  using  microlitho- 


graphy  processes  similar  to  those  used  to  make  computer  chips.  The  tips  protruding  from 


anisotropic  Si  etchant  is  then  used  to  free  the  cantilever  and  remove  the  portion  of  the  Si 


The  tips  of  the  cantilever  play  an  important  role  in  the  imaging  resolution  obtain- 


constant  while  the  cantilever  is  scanned  over  the  sample  surface.  The  topography  of  the  resonance  frequency  and  oscillation  amplitude  of  the  cantilever.  This  mode  operates  in 


material  property  information  for  a  sample. 


Table  2.1:  AFM  modes  used  to  measure  material  properties  such  as  surface  friction,  thermal 


depth,  and  contact  area,  to  the  material  properties  of  the  interacting  bodies,  such  as  Young* 


2.6.1  Hertzian  Contact  Model 

The  JKRS  model  is  an  extension  of  the  Hertzian  model,  which  incorporates  ad- 
Hertz  first  analyzed  the  elastic  contact  between  a  sphere  and  a  fiat  surface  in  1888  hesion  through  the  change  in  surface  energy  at  the  contact  of  the  two  spheres,  as  shown  in 

l']’I35l'  ffurtz  was  able  to  relate  the  contact  characteristics  of  the  applied  load,  indentation  2.9.  There  are  no  forces  between  the  spheres  when  they  are  not  in  contact  [17], [37].  This 


A  comparison  of  the  Maugis  contact  model  for  various  values  of  A  is  shown  in  Figure  2.11. 
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l.  Variations  in  the  conditions  allow  the  elastic  modulus  to  vary  within  the 


Figure  3.1  [40].  A  tip  with  a  small  radius  is  attached  to  the  end  of  the  cantilever. 


The  potential  energy  consists  of  two  parts:  one  due  to  bending  and  one  due  to  the  where  two  terms  are  retained  from  the  binomial  expansion.  Inserting  Equation  3.7  into 

axial  force  acting  on  the  beam.  Tire  potential  energy  density  due  to  the  gravitational  forces  Elation  3.6,  allows  the  combined  potential  energy  due  to  bending  and  the  axial  force  to 


ing  Equations  3,9,  3.11,  and  3.12  into  Equation  3.1  and  grouping  terms,  the 


assuming  harmonic  solutions  of  the  form,  w(xf  t)  —  ip(z)eMt,  which  allows  the  equation  of  By  equating  the  determinant  of  the  coefficients  to  zero,  the  characteristic  equation  governing 


C3(co6(,5L)  +  cosh(.dL))  +  C4(sm(/?I,)  +  &mh{j3L)) 


C^(cos\\(3L)  —  oos (3L))  -f  ^(siiih^L)  -  sin(5L)) 


is  assumed  that  the  cantilever  tip  and  the  sample  surface  deform  elastically  when  in  contact, 


lated  from  tip- sample  contact  models,  as  summarized  in  section  2.6.  The  research  for  this  The  derivation  of  the  equations  governing  the  torsional  vibrations  of  the  cantilevers 


where  the  torsional  wave  number,  7,  is  given  by 


The  stationary  of  Rayleigh’s  quotient  is  then  exploited  by  attempting  to  minimize  the 


moment  of  inertia,  J,  the  torsional  constant,  op,  and  the  area  polar  moment  of  inertia,  J, 
will  be  presented  for  each  model.  These  equations  will  then  be  used  with  the  Rayleigh-Ritz 


where  t  is  the  thickness  of  the  cantilever  and  W  is  the  width  of  the  cantilever.  For  torsional 
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the  sensitivity  for  the  first  mode  of  the  triangular  cantilever  is  slightly  greater  than  that  of 


equations  used  to  approximate  four  different  cantilever  shapes  were  presented  in  section 
4.3.  Example  modal  sensitivity  calculations  were  presented  to  demonstrate  the  usefulness  of 


X  10 _  .  I  this  selection  is  because  the  calculated  contact  resonance  frequencies  more  closely  match 


variations  in  the  thickness  of  the  cantilever. 


Table  5.3:  The  measured,  fcxp,  and  the  calculated,  featei  free  frequencies  for  the  rectangular 
NCS16  cantilever. 

mode  foxy  (kHz)  Scale.  (kHz) 

1 _ 178.0  175.5 


a  a 


e  a 

•&  p 

I  1 

I  I 

E  S 

4  1 


6 

1 
09 

l 
I 

=  o  .2 


eT 

a. 


II 

£ 

II 

sr 


£ 

p 

1 


p  1 

S  g 

3.  -5 

22  fl 


I  I 

a  ; 


p 

J8 


*3 

i 


?  ?  S 

8 
§ 


a, 

a 


ii 


II 

4 

2 


ci  .2 

1 


i  2 

o  js 


•S  JJ 


3 


6  » 


o 


!  '§ 

c§^  "8 


I  i 
s 


*3 

-a 


i 


"8 

.3 


13 

1 


O 

a) 

P 


n  »  U  V 

I  2  3  -o  3 


I 

as 

£ 


& 

z 


^  *s  o 


< 

S3 


1  » 

E  *3 

•2  2 


g 

a. 


I 

s 

p 

I 

g 


I 

•8 


1  1 


1 


P  P 


I  4 


s 

o 

B 


»  8 


8  *2 

2  S 

2  2 

o  =*. 

a  s 

1  8 

a  <n 

3  ii 


I 

•tt 

•8 


■2  *8 

8  | 


'S  aj 

J)  ** 

p  *6 


1 

1 


P 

i 


•i  ~ 

3  <2 


*0 

S 


I  i 


s 

I 


Table  5.5:  The  measured,  /exp,  and  the  calculated,  free  frequencies  for  the  NCSll  A 
cantilever. 


also  cause  the  calculated  contact  resonance  frequencies  to  differ  greatly  from  the  measured  5.3.2  Tbrsional  Free  Vibrations 


1063.0  1  046.4 


5,3.3  Flexural  Contact  Vibrations 


11127  M-fe 


xio 


*8  3  S 

.B  1  S 


o 

to 

3 


.o 

a 


i 


£ 


g  1 

I  I 

S  8 

1  I 


O  *-5  ~ 


£  3 


8  S 


ts 

I 

s 

a) 

*s 


S  ?  2  & 

£  C  CJ  'S3 


§  g 


s 

3 

I 

TJ 

3 

cC 

>> 

I 

TJ 


2 


1 


f— < 

£  i 


tis 


Jl 

I 


50c 


& 

*5 


1 


1 

8 


« 


it 

<S5C 


I  E 


S’ 


* 

*8 

§ 


I 


§ 

I 

a> 

P 


CO 

2 

53 


Jf 

Co. 

I 

| 


I 


i 

50. 


4 

50. 


-«JJ} 

5 

50. 


if 


7  -J? 

-i*  i 

"a*  50. 

a 

50.  -= 


—  ®  50. 

I  I  -!* 


4* 


i  fe 


&  1 


.2 

-a 


T3 
.2 
bfi 

.5  ftj 

f  ° 

i  I 


5  f 

1  i 
£  I 

S  1 


1 

i 

50. 


I 


I  1  I 


§  3 

5  as 


1  1  1 


I  .1 


s  «s 


I  £ 


■s'  °*  « 


53  I 


a 

*1  I 


4 

05 


1 

1 


3 

*a 


1 


I 


I  eg 


SO. 


4P 

50. 


Equation  5.8  is  used  to  calculate  stiffness  values  for  the  dLever  A  cantilever. 


tip-sample  interaction  is  calculated  from  Equation  2.2  as  E*  =  82.836  GPa.  The  radius  if  the 


eclmjque  the  error  jn  calculating  Es  is  around  20%  [43], [63].  This 


measuring  the  vibration  response  of  the  AFM  cantilevers  was  explained.  The  measured  free 


calculated  using  the  approximate  method,  compared  well  with  the  results  obtained  using 


the  elastic  modulus  of  a  sample  surface,  the  cantilever  used  for  the  measurements  must  be 


Force  Microscope  Cantilevers:  Theory  and  Experiment,”  Review  of  Scientific  Instru- 


AFAM  method,  surface  stiffnesses  of  these  samples  were  determined  and  some  images  of 

May,  2002 

the  stiffness  of  the  samples  were  taken,  the  results  of  which  are  discussed  and  analyzed. 


3.8  SEM  image  of  silicon  calibration  grating  TGT01  of  an  array  of  sharp  tips  4.13  Frequency  shift  due  to  different  applied  loads  for  cantilever  NSC-16  on  sample 

characterized  by  strict  symmetry  of  tip  sides,  small  cone  angles  and  small  Si<lll>,  voltage  amplitude  IV . 
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ultrasonic  excitations  -.'induced  in  the  sample  or  base  of  the  cantilever.  In  this  case, 
the  contact  between  the  cantilever  tip  and  the  specimen  causes  the  natural  frequencies  of 


effects.  AFMs  can  achieve  a  resolution  of  up  to  1  mn,  and  unlike  electron  microscopes, 


fields  may  also  be  applied  intentionally  between  the  tip  and  the  sample. 


Figure  2.5:  AFM  3D  topography  image  of  a  CrBN  thin  film.  atoms  reaches  a11  equilibrium  distance,  which  is  about  a  couple  of  angstroms.  If  the  atoms 


topography  by  keeping  the  cantilever  deflection  constant.  In  this  case,  the  image  is  generated 
from  the  scanner’s  motion.  With  the  cantilever  deflection  held  constant,  the  total  force 
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■  ^  reality,  the  cross  section  of  most  cantilevers  can  deviate  from  a  perfect  rectangular  I  The  deflection  i/j(x,  i)  is  defined  as  the  cantilever  position  relative  to  its  initial 


Figure  2.7:  A  schematic  of  the  cantilever  in  contact  with  a  sample. 


High-resolution  acoustic  imaging  can  be  a  powerful  tool  for  materials  research. 


with  increased  excitation  amplitude  or  decreased  static  cantilever  load  the  contact 


environment  (normally  on  a  silicon  chip).  Microfabrication  technology  enables  fabrication 


devices  can  accomplish  complicated  functions.  Current  MEMS  applications  include,  but 


In  this  chapter,  the  experimental  setup  developed  is  discussed  in  detail. 
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for  most  cantilevers,  the  force  constant  is  less  than  1  N/m.  Cantilevers  may  be  divided  into 
two  categories,  either  soft  or  hard,  in  terms  of  their  force  oonstant.  Soft  cantilevers  which 
have  a  force  constant  less  than  1  N/in  are  usually  used  for  experiments  on  soft  samples 


After  mounting  the  desired  cantilever  on  the  probe  cartridge,  the  head  of  the  AFM 


vertical  and  lateral  calibration  of  the  AFM  scanner,  detect  vertical  and  lateral  nonlinearity 


cantilever),  the  program  outputs  are  presented  in  Appendix  A. 


the  tip  is  very  sharp  and  new,  taking  a  high  resolution  SEM  image  for  determining  tip  radius 


Figure  3.23:  An  SEM  image  the  tip  of  a  d-lever  B  cantilever. 


3.4  Surface  stiffness  calculation 

Table  3.3:  Some  mechanical  properties  of  the  tip  and  the  sample. 

Using  the  Hertzian  contact  theory  described  earlier  in  Section  2.2.4,  it  is  possible  —  0,180  Pois8Qn  coefficient  for  the  sample  I 
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surface  stiffness  map  of  the  image  is  presented.  In  the  last  section,  some  nonlinear  effects 


AFAM  images  of  the  surface  of  poiysilioon  MEMS  samples  are  presented,  as  is  the  related 


vibrating  cantilever.  The  damping  can  be  measured  using  the  peak  width,  where  a  wider 


4.1.2  d- lever  B  cantilever 


to  determine  surface  stiffness. 


Cantilever  NSC-10  frequency  shift  response  Fbr  the  NSC-16  cantilever  the 
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Figure  4.12:  d-lever  B,  repeatability  of  the  data  for  applied  force  120  nN.  SiSn^TOlt^^litude  IV  *°  diffetent  '0adS  f°r  cantilever  NSC-16  °»  ^"iple 


with  increasing  applied  forces  to  the  cantilever.  According  to  Figure  4.16,  the  first  natural 


from  the  lock-in  amplifier.  When  the  cantilever 


taken  using  the  first  mode  /=  837  kHz. 


«  c  "S  s 

**  N  W  ^ 


5  a  a  * 

00  w  £ 

CM  CO  -C 


Z  Z  Z  Z 

c  c  c  c 

o  o  o  o 

■>tf-  -'a-  ■<»- 


z  z  z  z 

c  s  t=  c 
o  o  o  o 
(O  ffl  o  o 


IUUJUlUJUJUIUJtUUl|^ 

OOOOOOOOOq 

ooco^rcNjpqoooQ 

T-"  T-’  r-  r-  T-  CO  O  T  N  g 


LULUUJUiUJLUlUlUUJUJyj 

OOOOOOOOOOq 

oa>oo^rcNippppp0 

oi  r  r  ^  r  r  fl)  c  V  (N  o 


Figure  4.28:  Nonlinear  repeatability  on  d-lever  B,  applied  force  60  nN,  applied  voltage  to 
the  transducer  5V. 
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Figure  4.30:  Nonlinear  repeatability  on  d-Iever  B,  applied  force  100  nN,  applied  voltage  to 
the  transducer  5V. 


sensitivities  using  other  cantilevers  seems  a  logical  next  step.  Specially,  in  order  to  find 


[7]  L.  Meirovitch,  Analytical  Methods  in  Vibrations  (The  Macmillan  Company,  1975). 
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